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Abstract 

Within a light-cone quantum-chromodynamics dipole formalism based on the Green 
function technique, we study nuclear shadowing in deep-inelastic scattering at small 
Bjorken xsj 0.01. Such a formalism incorporates naturally color transparency and 
coherence length effects. Calculations of the nuclear shadowing for the qq Fock compo- 
nent of the photon are based on an exact numerical solution of the evolution equation 
for the Green function, using a realistic form of the dipole cross section and nuclear den- 
sity function. Such an exact numerical solution is unavoidable for xbj ^ 10~^, when a 
variation of the transverse size of the qq Fock component must be taken into account. 
The eikonal approximation, used so far in most other models, can be applied only at 
high energies, when xbj ^ 10~^ and the transverse size of the qq Fock component is 
"frozen" during propagation through the nuclear matter. At xbj < 0.01 we find quite a 
large contribution of gluon suppression to nuclear shadowing, as a shadowing correction 
for the higher Fock states containing gluons. Numerical results for nuclear shadowing 
are compared with the available data from the E665 and NMC collaborations. Nuclear 
shadowing is also predicted at very small XBj corresponding to LHC kinematical range. 
Finally the model predictions are compared and discussed with the results obtained from 
other models. 
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1 Introduction 



Nuclear shadowing in deep-inelastic scattering (DIS) off nuclei is usually studied via nuclear 
structure functions. In the shadowing region of small Bjorken Xbj ^ 0.01 the structure function 
F2 per nucleon turns out to be smaller in nuclei than in a free nucleon (see the review [1], for 
example). This affects then the corresponding study of nuclear effects, mainly in connection 
with the interpretation of the results coming from hadron-nucleus and heavy ion experiments. 

Nuclear shadowing, intensively investigated during the last two decades, can be treated 
differently depending on the reference frame. In the infinite momentum frame of the nucleus it 
can be interpreted as a result of parton fusion [21 [3l IH [5], leading to a reduction of the parton 
density at low Bjorken xbj- In the rest frame of the nucleus, however, this phenomenon looks 
like nuclear shadowing of the hadronic fluctuations of the virtual photon, and occurs due to 
their multiple scattering inside the target P[71[8l[9l[T0l[IIl[T2l[T3l[Tll[T5l[T6l[I^ Although 
these two physical interpretations are complementary, the one based on the rest frame of the 
nucleus is more intuitive and straightforward. 

The dynamics of nuclear shadowing in DIS is controlled by the effect of quantum coherence, 
which results from the destructive interference of amplitudes for which the interaction takes 
place on different bound nucleons. Taking into account the \qq) Fock component of the photon, 
quantum coherence can be characterized by the lifetime of the qq fluctuation, which in turn 
can be estimated by relying on the uncertainty principle and Lorentz time dilation as, 

where u is the photon energy, is photon virtuality and Mgg is the effective mass of the qq pair. 
This is usually called coherence time, but we also will use the term coherence length (CL), since 
light-cone kinematics is assumed, Ic = t^. The CL is related to the longitudinal momentum 
transfer by qc = l/lc- Notice that for higher Fock states containing gluons \qqG), \qq2G), ... , 
the corresponding effective masses are larger than Mgg. Consequently, these fluctuations have a 
shorter coherence time than the lowest \qq) state. The effect of CL is naturally incorporated in 
the Green function formalism, which has been already applied to DIS, Drell-Yan pair production 
[121 [TTl [IS] , and vector meson production [201 EI] (see also the next Section). 

In the present paper nuclear shadowing in DIS will be treated using again the Green function 
approach. Such a quantum mechanical treatment requires to solve the evolution equation for 
the Green function. Usually, for simplicity this equation is set up in a such way as to obtain 
the Green function in an analytical form (see [121 [I?], for example), which requires, however, 
to implement several approximations into a rigorous quantum-mechanical approach, like a 
constant nuclear density function fl36|) and a specific quadratic form fl35|) of the dipole cross 
section. The solution obtained in a such way is the harmonic oscillator Green function [22] 
(see also Eq. f|T6|) ). usually used for calculation of nuclear shadowing [HI HTl [23]. Then the 
question about the accuracy of the predictions for nuclear shadowing using such approximations 
naturally arises. 

In the process of searching for the corresponding answer, in 2003 the evolution equation for 
the Green function was solved numerically for the first time in ref. [H] . This allowed to exclude 
any additional assumptions and avoid supplementary approximations, which caused theoretical 
uncertainties. The corresponding predictions for nuclear shadowing in DIS at small Xsj, based 
on the exact numerical solution of the evolution equation for the Green function [IB], showed 
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quite a large difference in comparison with approximate calculations [TSl [IZ] obtained within 
the harmonic oscillator Green function approach, in the kinematic region when Ic < Ra [Ra is 
the nuclear radius). However, no comparison with data was performed using this path integral 
technique based on an exact numerical solution of the two-dimensional Schrodinger equation 
for the Green function. This is one of the main goals of the present paper. Such a comparison 
with data provides a better baseline for future studies of the QCD dynamics, not only in DIS 
off nuclei but also in further processes occurring in lepton (proton)-nucleus interactions and in 
heavy-ion collisions. 

The calculations of nuclear shadowing in DIS off nuclei presented so far within the light- 
cone (LC) Green function approach [191 [13 HE] were performed assuming only qq fluctuations 
of the photon, and neglecting higher Fock components containing gluons and sea quarks. The 
effects of higher Fock states are included in the energy dependence of the dipole cross section, 
-sjll- However, as soon as nuclear effects are considered, these Fock states \qqG), \qq2G) 

lead to gluon shadowing (GS), which for simplicity has been neglected so far when the model 
predictions were compared with experimental data. The contribution of the gluon suppression 
to nuclear shadowing represents a shadowing correction for the multigluon higher Fock states. 
It was shown in ref. [24j that GS becomes effective at small xbj < 0.01. The present available 
experimental data cover the shadowing region ~ 0.0001 < XBj ^ 0.01, and therefore the 
contribution of GS to the overall nuclear shadowing should be included. This is a further goal 
of the present paper. 

Different (but equivalent) descriptions of GS are known, depending on the reference frame. 
In the infinite momentum frame of the nucleus it looks like fusion of gluons, which overlap in the 
longitudinal direction at small xsji leading to a reduction of the gluon density. In the rest frame 
of the nucleus the same phenomenon looks as a specific part of Gribov's inelastic corrections 
[25] . The lowest order inelastic correction related to diffractive dissociation j* N X N [26] 
contains PPR and PPP contributions (in terms of the triple- Regge phenomenology, see [?F]). 
The former is related to quark shadowing, while the latter, the triple- Pomeron term, corresponds 
to gluon shadowing. Indeed, only diffractive gluon radiation can provide the Mx dependence 
do'dd/dMx oc 1/Mx of the diffractive dissociation cross section. In terms of the light-cone 
QCD approach the same process is related to the inclusion of higher Fock components, \qqnG), 
containing gluons [28]. Such fluctuations might be quite heavy compared to the simplest \qq) 
fluctuation, and therefore have a shorter lifetime (see Eq. ([T])), and need higher energies to be 
relevant. 

Calculations of the GS contribution to nuclear suppression have been already performed 
within the light-cone QCD approach, for both coherent and incoherent production of vector 
mesons [2Ul UJ] . and also for production of Drell-Yan pairs [T7]. They showed (except for 
the specific case of incoherent production of vector mesons) that GS is a non-negligible effect, 
especially for heavy nuclear targets at small and medium values of photon virtualities < a 
few GeV^ and at large photon energies u. This is another reason to include the effect of GS 
for the calculation of nuclear shadowing, especially for making more realistic comparison of the 
predictions with experimental data. 

Notice also that by investigating shadowing in the region of small xbj < 0.01 we can safely 
omit the nuclear antishadowing effect assumed to be beyond the shadowing dynamics [H |9] . 

The paper is organized as follows. In the next Section [2] we present a short description of the 

^Here r represents the transverse separation of the qq photon fluctuation and s is the center of mass energy 
squared. 
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light-cone dipole phenomenology for nuclear shadowing in DIS, together with the Green function 
formalism. In Section [3] we discuss how gluon shadowing modifies the total photoabsorption 
cross section on a nucleus. In Section H] numerical results are presented and compared with 
experimental data, and also with the results from other models, in a broad range of xbj- Finally, 
in Section [5] we summarize our main results and discuss the possibility of future experimental 
evidence of the GS contribution to the overall nuclear shadowing in DIS at small values of xbj- 



2 Light-cone dipole approach to nuclear shadowing 

In the rest frame of the nucleus the nuclear shadowing in the total virtual photoabsorption 
cross section dtot^ixsj, Q"^) (or in the structure function F^{xBj, Q'^)) can be decomposed over 
different Fock components of the virtual photon. Then the total photoabsorption cross section 
on a nucleus can be formally represented in the form 

c^tot^ixBj, Q^) = A a1!*f{xBj, Q^) - Aatot{xBj, Q^) , (2) 

where 

AatotixBj, Q^) = Aatotiqq) + Acxtot{qqG) + Aatotiqq'^G) + • . (3) 
Here the Bjorken variable xbj is given by 

''''' 2mr,u^Q^ + s' 

where s is the 7*-nucleon center of mass (cm.) energy squared, m^r is mass of the nucleon, 
and CTtot^ {x Bj , Q"^) in is total photoabsorption cross section on a nucleon 
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(^M^{xBj,Q^) = J d^r da '^qq{f,a,Q'^) aqq{f,s). (5) 

In this last expression <Jqq{r, s) is the dipole cross section, which depends on the qq transverse 
separation r and the cm. energy squared s, and "^gqir, a, Q'^) is the LC wave function of the qq 
Fock component of the photon, which depends also on the photon virtuality and the relative 
share a of the photon momentum carried by the quark. Notice that XBj is related to the cm. 
energy squared s via Eq. (jlj). Consequently, hereafter we will write the energy dependence of 
variables in subsequent formulas also via an x^j-dependence whenever convenient. 

The total photoabsorption cross section on a nucleon target ^ contains two ingredients. 
The first ingredient is given by the dipole cross section (Jqq{r, s), representing the interaction of 
a qq dipole of transverse separation r with a nucleon [29j. It is a flavor independent universal 
function of r and energy, and allows to describe various high energy processes in an uniform way. 
It is also known to vanish quadratically (Tqq{r, s) oc as r — >• 0, due to color screening (property 
of color transparency [2H1 EOl EI]), and cannot be predicted reliably because of poorly known 
higher order perturbative QCD (pQCD) corrections and nonperturbative effects. However, 
it can be extracted from experimental data on DIS and structure functions using reasonable 
parametrizations, and in this case pQCD corrections and nonperturbative effects are naturally 
included in aqqir^s). 

There are two popular parameterizations of aqqif^s): GBW presented in [32], and KST 
proposed in [21]. Detailed discussions and comparison of these two parametrizations can be 
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found in refs. [221 EHl [IB]- Whereas the GBW parametrization cannot be apphed in the non- 
perturbative region of Q^, the KST parametrization gives a good description of the transition 
down to the hmit of real photoproduction, = 0. Because we will study the shadowing region 
of small XBj ^ 0.01, where available experimental data from the E665 and NMC collaborations 
cover small and moderate values of < 2 ^ 3 GeV^, we will prefer the latter parametrization. 

The KST parametrization [2l] has the following form, which contains an explicit dependence 
on energy. 




(6) 



exp 

The explicit energy dependence in the parameter ao{s) is introduced in a such way that it 
guarantees that the correct hadronic cross sections is reproduced, 

Ms) = ^im fi + |#^l , (7) 



{r. 
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where crtJl{s) = 23.6 (s/so)° °'^^ + 1.432 {s/ Sq)~^''^^ mh, which contains the Pomeron and Reggeon 
parts of the up total cross section ^33j, and Ro{s) = 0.88 (s/sq)"'^^^ fm, with A = 0.28 and where 
So = 1000 GeV^ is the energy-dependent radius. In Eq. ([?[) {rli^)n = 0.44 fm^ is the mean pion 
charge radius squared. The form of Eq. ([6|) successfully describes the data for DIS at small x Bj 
only up to ~ 10 GeV^. Nevertheless, this interval of is sufficient for the purpose of the 
present paper, which is focused on the study of nuclear shadowing at small xsj ^ 0.01 in the 
kinematic range 

Q2 < 4 GeV^ covered by available E665 and NMC data. 
However, as we will present the predictions for nuclear shadowing at very small xbj down to 
10~^ accesible by the prepared experiments at LHC and at larger values of > 10 GeV^, we 
will use also the second GBW parametrization [32j of the dipole cross section. 

The second ingredient of cr^^t^ {x bj , Q"^) in ([Sj) is the perturbative distribution amplitude 
( "wave function" ) of the qq Fock component of the photon. For transversely (T) and longitu- 
dinally (L) polarized photons it has the form [Ml [3S1 fTU] : 



^^^\r,a,Q') = ^^^^ Z^xO^'^^xMer) (8) 



where x and x the spinors of the quark and antiquark respectively, Zg is the quark charge, 
Nc = 3 is the number of colors, and Ko{er) is a modified Bessel function with 

= a{l~a)Q^ + ml , (9) 
where ruq is the quark mass. The operators O^'^ read, 

O'^ = nig a ■ e + i{l - 2a) {a ■ n) (e ■ Vr-) + (a x e) ■ , (10) 

= 2Qa{l- a){a -n) . (11) 

Here acts on the transverse coordinate r, e is the polarization vector of the photon, n is a unit 
vector parallel to the photon momentum, and a is the three vector of the Pauli spin-matrices. 

The distribution amplitude Eq. ([8]) controls the transverse qq separation with the mean 
value 

{r)r^- = — . 12 
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For very asymmetric qq pairs with a or (1 — a) < w^jQ^ the mean transverse separation 
(r) ~ l/wiq becomes huge, since one must use current quark masses within pQCD. A pop- 
ular recipe to fix this problem is to introduce an effective quark mass me// ~ ^qcd-, which 
represents the nonperturbative interaction effects between the q and q. It is more consistent 
and straightforward, however, to introduce this interaction explicitly through a phenomenology 
based on the light-cone Green function approach, and which has been developed in YIa 



The Green function G'gg(r'^, ^2; -^i) describes the propagation of an interacting qq pair 
between points with longitudinal coordinates Z\ and Z2 and with initial and final separations 
f\ and r2. This Green function satisfies the two-dimensional Schrodinger equation, 



dZ2 



2 u a (1 — a] 



+ Vqq{Z2,f2,a) 



Gqq{f*2,Z2;fl,Zi) 



(13) 



with the boundary condition 



Gqq{r2i Z2', ri, Zi] 



(14) 



In Eq. f fT3|) u is the photon energy and the Laplacian A,, acts on the coordinate r. 

We start with the propagation of a qq pair in vacuum. The LC potential Vqq{z2, r^, a) in f|T3l) 
contains only the real part, which is responsible for the interaction between the q and q. For the 
sake of simplicity we use an oscillator form of this potential. Although more realistic models 
for the real part of the potential are available [36l [37] , however, solution of the corresponding 
Schrodinger equation for the light-cone Green function is a challenge. Analytic solution has 
been known so far only for the oscillator potential. Otherwise one has to solve the Schrodinger 
equation numerically, which needs a dedicated study. 

On the other hand, important is the mean qq transverse separation which is fitted to diffrac- 
tion data. Any form of the potential must comply with this condition. The same restriction 
is imposed on the quark-gluon Fock states. The mean quark-gluon separation, which matters 
for shadowing, is fixed by high-mass diffraction data and should not be much affected by the 
choice of a model for the potential. 



ReVqq{z2,f^2,a) 



a'^{a) 7^2 ^ 
2 u a(l — a) 



(15) 



one can solve then two-dimensional Schrodinger equation (fT3|l analytically, and the solution is 
given by the harmonic oscillator Green function [38] 



Gqq{r 2, Z2',ri, Z\] 



2 Ti i siniuj Az" 



exp 



sinftu Az) 



[r\ + Tg) cos(co' Az) — 2 fl ■ f*2 



xexp 



ie^Az 



2 u a (1 ~ a) 



(16) 



where Az = Z2 — Zi, and 



0^(0;) 



u a{l — a) 

The shape of the function a{a) in Eq. f|T5l) will be discussed below. 



(17) 
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The probabihty amphtude to find the qq fiuctuation of a photon at the point Z2 with 
separation r, is given by an integral over the point Zi where the qq is created by the photon 
with initial separation zero, 



An E ail - a) 



J dZi (x O^'^X) Ggq{f, Z2]fl,Zi) 



ri=0 



:i8) 



The operators O^'^ are defined by Eqs. (JTOl) and (fTTl) . and here they act on the coordinate rl. 
If we write the transverse part as 



X O'^X = X^c ^ ■ ex + X - '^a) {a ■ n) e + {a ^ e)] X - ^r = E + F - V r 
then the distribution functions read, 

^Ur,a) = Zqy/^ [E$o(e,r,A) +F$i(e,r,A) 



where 



2 Zq^/a^Q a{l-a)x^ -nx '^'o(e, r. A) 

_ 2a2(a) 
A — — . 



The functions $0,1 iii Eqs. fl20l) and fl2Tl) are defined as 



1 °° 

^o(e,r,A) = ^Jdt 



A 



sh(At) 



exp 



Xe^r 



2^2 



$i(e,r,A) = — /rft 
Stt 



00 


A 


2 


/* 




exp 


_sh(At)_ 



4 

AeV2 
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cth(At) - 1 
cth(At) - 1 



(19) 

(20) 
(21) 

(22) 

(23) 
(24) 



where sh{x) and cth{x) are the hyperbolic sine and hyperbolic cotangent, respectively. 

Note that the q — q interaction enters in Eqs. (|20l) and (!2T1) via the parameter A defined in 
Eq. ( I22I) . In the limit of vanishing interaction A — (i.e. — > 00, a is fixed, a 7^ or 1) 
Eqs. ( I20I) - (1211) produce the perturbative expressions of Eq. ([H]). 

With the choice a^(a;) oc a{l — a), the end-point behavior of the mean square interquark 
separation is (r^) oc l/a(l — a), which contradicts the idea of confinement. Following [23] we 
fix this problem via a simple modification of the LC potential, 

a^{a) = al + 4ala{l - a) . (25) 

The parameters qq and oi were adjusted in [23] to data on total photoabsorption cross section 
[5^ 1^ . diffractive photon dissociation, and shadowing in nuclear photoabsorption reaction. 
The results of our calculations vary within only 1% when oq and ai satisfy the relation. 



1 



1-1^ (0.112)2 (^gy2 
1.15 /n 1 (ir:\2 



(0.165)2 GeV" , 



(26) 



where v takes any value < f < 1. In view of this insensitivity of the observables we fix the 
parameters at f = 1/2. We checked that this choice does not affect our results beyond a few 
percent uncertainty. 
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The matrix element contains the LC wave function squared, which has the following 
form for T and L polarizations, in the limit of vanishing interaction between q and q, 



2Nr:a 



Nf 



C Ltem 



and 



(2vr) 



E 



ZfQ a [l-a) Ko{er) , 



where Ki is the modified Bessel function. 



d 



Ki{z) = -—Ko{z) 



(27) 



(2^ 



(29) 



If one includes the nonperturbative q — q interaction, the perturbative expressions fl27j) and 
( !28l) should be replaced by: 



Nf 



Cpi(f,a,Q2) 



2 Nc aem E [m} ^^(e, r. A) + [a^ + (1 - a)^] |$i(e, r. A) 



and 



AT, 



(30) 



(31) 



/=i 

Notice that in the LC formalism the photon wave function contains also higher Fock states 
\qq), \qqG), \qq2G), etc., but its effect can be implicitly incorporated into the energy dependence 
of the dipole cross section o"gg(r, s), as is given in Eq. ([5]). The energy dependence of the dipole 
cross section is naturally included in the realistic KST parametrization Eq. ([6]). 

Now we will continue with our discussion of DIS on nuclear targets, and will study the 
propagation of a qq pair in nuclear matter. Some work has already been done in this direction. 
In fact, the derivation of the formula for nuclear shadowing, keeping only the first shadowing 
term in Eq. ([2]), AatotixBjjQ"^) = 'Aortot{<l<l) , can be found in [H]. This term represents the 
shadowing correction for the lowest qq Fock state, and has the following form 



Aatot{xBj,Q'^ 



Re / d% 



dzi pA{b,zi) / dz2 pA{b,Z2) / da A{zi,Z2,a) , (32) 



with 

A{zi,Z2,a) = J d^r2 ^^^(r^, a, Q^) agq{r2,s) J rfVi Ggq{f'2, Z2;fl, zi) agg{ri,s) ^^^(rl, a, Q^) . 

(33) 

When nonpertubative interaction effects between the q and q are explicitly included, one should 
replace in Eq. <il gg{f, a , Q^) =^ ^„pt(f,a,Q2) and <i!t^{r,a,Q^) =^ ^^^.(f, a, Q^). 

In Eq. ([32D PAip,z) represents the nuclear density function defined at the point with longi- 
tudinal coordinate z and impact parameter h. 

The shadowing term /S.atot{xBj,Q'^) = Aatotiqi) in ([2]) is illustrated in Fig. [TJ At the point 
Zi the initial photon diffractively produces the qq pair (7*A^ — >■ qqN) with transverse separation 
fi. The qq pair then propagates through the nucleus along arbitrary curved trajectories, which 
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Figure 1: A cartoon for the shadowing term AatotixBjyQ"^) = Ao'toti'i'i) 
in ([2]). Propagation of the qq pair through the nucleus is described by 
the Green function Ggq(r'2, ^2; ^, -Zi), which results from the summation 
over different paths of the qq pair. 



are summed over, and arrives at the point Z2 with transverse separation r^. The initial and 
final separations are controlled by the LC wave function of the qq Fock component of the 
photon '^qq{r, a, Q"^). During propagation through the nucleus the qq pair interacts with bound 
nucleons via the dipole cross section crgg(r, s), which depends on the local transverse separation 
r. The Green function Gqq{f^, Z2;fl,Zi) describes the propagation of the qq pair from zi to Z2- 
Describing the propagation of the qq pair in a nuclear medium, the Green function 
Gqq{f2, Z2]fi,zi) satisfics again the time-dependent two-dimensional Schrodinger equation f|T3|) . 
However, the potential in this case acquires in addition an imaginary part. This imaginary part 
of the LC potential Vqq{z2, r2, in Eq- dlSD is responsible for the attenuation of the qq photon 
fluctuation in the medium, and has the following form 

ImVqq{z2, r-, a) = p^{b, Z2) . (34) 

As was already mentioned above, the analytical solution of Eq. ( |T3i) is known only for the 
harmonic oscillator potential Vqq{r) oc r^. Consequently, in order to keep such an analytical 
solution one should also use a quadratic approximation for the imaginary part of Vgq{z2, r2, «), 
i.e. 

aqqir,s) = C{s)r\ (35) 

and uniform nuclear density 

PA{b,z)=poe{R\-b'-z^). (36) 

In this case the solution of Eq. ( JT3l) has the same form as Eq. ( JT6il . except that one should 
replace uj =^ Q and 0^(0) =^ b{a), where 

^ _ b{a) _ \l a^(a) - i PA^b, z)ua{l-a) C{s) ^^^^ 
z/a(l — a) 



v a 



[I -a) 



9 



The determination of the energy dependent factor C{s) in Eq. fl35p and the mean nuclear 
density po in Eq. fl36l) can be reahzed by the procedure described in [T71 |23l [18] , and will be 
discussed below. 

Investigating nuclear shadowing in DIS one can distinguish between two regimes, depending 
on the value of the coherence length: 

(i) We start with the general case when there are no restrictions for Ic- If Ic ~ Ra one has to 
take into account the variation of the transverse size r during propagation of the qq pair through 
the nucleus, which is naturally included using a correct quantum-mechanical treatment based on 
the Green function formalism presented above. The overall total photoabsorption cross section 
on a nucleus is given as a sum over T and L polarizations, a"'*'^ = cxj. ^ + e' assuming that 
the photon polarization e' = 1. If one takes into account only the qq Fock component of the 
photon, the full expression after summation over all flavors, colors, helicities and spin states 
becomes HJl 



+ 



C «em 



J r poo poo r r 

^ i?e / d^b / dzi / dz2 da d\i / d' 

r-^ J J ~00 J Zl Jo J J 



(2^) f 

X pA{b, Zl) pA{b, Z2) (Tgq{r2, s) 0-55 (ri, s 



r2 



X 



a 



+ 



Ti r2 



(38) 



mj + AQ' a' {I - ay KQ^eri) KQ{er2)\ Gqq{r2. Z2]fi, z^) 



Here '^'^^(f, a, Q"^) are the absolute squares of the LC wave functions for the qq fluctuation 
of T and L polarized photons, summed over all flavors, and with the form given by Eqs. (1271) 
and (!28l) . respectively. 

If one takes into account the nonperturbative interaction effects between q and q of the 
virtual photon the expression for a"'*^{xBj,Q'^) Eq. (l38l) takes the following form: 



cr, 



7* A 
npt 



{xBj,Q^ 



A Crlpt{xBj, Q'^) - A CTnptixBj, Q'^ 



A / 



da (Jqq{r, s) 



+ 



Nf 



Nc 0!em Zj Re / d^b / dzi / dz2 da d^ri I d^r2 

J— I J —00 J Zl Jo J J 

X pA{b, Zl) pA{b, Z2) (Tqq{r2, s) (Jqq{ri, s) 



X 



a^ + {l-a)^ $i(e,ri,A) ■ <l>i(e,r2,A) 



(39) 



+ 



mj + 4 Q2 _ ^)2 j ^ A) $o(e , ra, A)| Gg-q(r"^, Z2] r"l, Zi) 



where 



'^nyt (^) ct) Q^) are now given by Eqs. (l30l) and (l3Ti) . respectively. 
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(ii) The CL is much larger than the mean nucleon spacing in a nucleus {Ic ^ Ra), which is 
the high energy limit. Correspondingly, the transverse separation r between q and q does not 
vary during propagation through the nucleus (Lorentz time dilation). In this case the eikonal 
formula for the total photoabsorption cross section on a nucleus can be obtained as a limiting 
case of the Green function formalism. Indeed, in the high energy limit z/ — » oo, the kinetic term 
in Eq. ( fT3l) can be neglected and the Green function reads 



Gqq{b; rl, Z2; n, 2i)|i.^oo = S{r2 - n) exp 



2 ^991^2, S J 



22 



dz pA{b,z) 



21 



(40) 



Including nonperturbative interaction effects between q and q, after substitution of the expres- 
sion (HOll into Eq. (l39l) . one arrives at the following results: 



cr. 



7* A 
npt 



{xBj,Q^) = 2 J d% J d\ £ da(^l-exp 



-^^ggir, s)TA{b) 



x2Nr:a 




a 



$i(e,r,A) 



(41) 



+ 



mj + AQ'a'il-ay $^(e,r,A) 



where 



TAib) 



dz pA{b,z) 



(42) 



is the nuclear thickness calculated with the realistic Wood-Saxon form of the nuclear density, 
with parameters taken from [43]. 

At the photon polarization parameter e' = 1 the structure function ratio / is related 
to nuclear shadowing R{A/N and can be expressed via a ratio of the total photoabsorption 
cross sections 



Fi'ixB^Q' 



AR{A/N) 



7*A/ 



(43) 



where the numerator on the right-hand side (r.h.s.) is given by Eq. (!39|) . whereas the denomi- 
nator can be expressed as the first term of Eq. (139!) divided by the mass number A. 

As we already mentioned above, an explicit analytical expression for the Green function 



Gqq{f2, Z2',ri,zi) (fT6l) cau be found only for the quadratic form of the dipole cross section ( |35l) . 
and for uniform nuclear density function ( |36l) . It was shown in refs. p!9l [23l [T71 [18] that such 
an approximation gives results of reasonable accuracy, especially at small xbj ^ 10"'' and 
for heavy nuclei. Nevertheless, it can be even more precise if one considers the fact that the 
expression fHTl) in the high energy limit can be easily calculated using realistic parametrizations 
of the dipole cross section (see Eq. for the KST parametrization and ref. [32] for the GBW 
parametrization) and a realistic nuclear density function pA^b, z) [13]. Consequently, one needs 
to know the full Green function only in the transition region from non- shadowing {xBj ~ 0.1) 
to a fully developed shadowing given when coherence length Ic ^ -Ra, which corresponds to 
XBj ^ 10""^ depending on the value of Q^. Therefore, the value of the energy dependent factor 
C(s) in Eq. ( !35l) can be determined by the procedure described in refs. [171 [23], [20] . According 
to this procedure, the factor C(s) is adjusted by demanding that calculations employing the 
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approximation (153]) reproduce correctly the results for nuclear shadowing in DIS based on the 
realistic parametrizations of the dipole cross section Eq. ([6]) in the limit Ic ^ Ra, when the 
Green function takes the simple form (HOj) . Consequently, the factor C{s) is fixed by the relation 



JdH Jd'^r 




,Q2)p {i-exp 


'-lC{s)r'TA{b) 


} 


J d'^ r 







JdH Jd'^r 




2 r 

1 1 — exp 


-\(^qq{r,s)TA{h) 


} 






2 



(44) 



Correspondingly, the value po of the uniform nuclear density (l36l) is fixed in an analogous way 
using the following relation 



dH 



exp 



cTo poJR\- 62 



1 - exp l--aoTA{b) 



(45) 



where the value of po was found to be practically independent of the cross section ctq, when this 
changed from 1 to 50 mb [23]. Such a procedure for the determination of the factors C{s) 
and Po was applied also in refs. [20l [21], in the case of incoherent and coherent production of 
vector mesons off nuclei. 

In order to remove the above mentioned uncertainties the evolution equation for the Green 
function was solved numerically for the first time in ref. [TH]. Such an exact solution can be 
performed for arbitrary parametrization of the dipole cross section and for realistic nuclear 
density functions, although the nice analytical form for the Green function is lost in this case. 

In the process of numerical solution of the Schrodinger equation ([T3l) for the Green function 
Gqq{f2, Z2;'fi, Zi) wlth the initial condition ([T^ . it is much more convenient to use the following 
substitutions [18] 



go{r^2, Z2] zi, X) = d^ri <l>o(e , ri. A) agq(ri, s) ^^^(r^, Z2] rl, Zi) , 



and 



r2 



gi{r^2, Z2; zi, X) = / d'^ri$i{e,ri,X)aqg{ri,s)Gqq{r2,Z2;ri,Zi) . 



(46) 



(47) 



After some algebra with Eq. (IT^ these new functions go{f'2, Z2] Zi, X) and (?i(r^, Z2] ^i. A) can 
be shown to satisfy the following evolution equations 



9o[r2^Z2]Zi,X) 

dZ2 



and 



9iir2, Z2; zi, X) 

dz2 



2 Hqq 



1 



1 d 



52 



dr\ r2 dr2 



+ Vqq{z2, fa, a) \ go{f*2, ^2; ^1, A) (48) 



^ 2 flqq 

with the boundary conditions 

9oir^2, Z2; zi, X)\ 



92 Id I' 



dr2 r2 dr2 



+ Vqq{z2,r2,a) > fi-i (r^ , ^2 ; 2^1 , A) 



$o(e,r2,A) aqq{r2,s) 



(49) 
(50) 
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and 

9iir^2, Z2; zu X)\^^=^^ = <li(e,r2, A)crgq(r2,s) , 
where $i(e , r, A) is connected with $i(e , r, A) by the following relation: 

r ~ 



$i(e,r, A) = -<l>i(e,r,A) 



In Eqs. (SH]) and dH]) the quantity 



Hqq = u a{l — a) 



plays the role of the reduced mass of the qq pair. 

Now the expression (l39l) for total photoabsorption cross section on a nucleus reads 



■y*N 



A I d r J da (Jqq{r, s) 

Nf 



+ 



f r poc poo pi p ^ 

— 3 ^ ^/ Re / d^b / dzi / dz2 / da d' 

j^—l J— 00 J z\ Jo J 

X pA{b, Zi) pA{b, Z2) (Tqq{r2, s) 

+ (1 - af $i(e,r2, A)5fi(r^,Z2;^i, A) 



r2 



X 



(51) 



(52) 



(53) 



(54) 



+ 



mj + 4 Q2 ^2 _ ^^^^ ^ ^-^ ^ 1 _ 



Notice that this equation explicitly includes nonperturbative interaction effects between q and 
q. Details of the algorithm for the numerical solution of Eqs. (I48l) and (l49l) can be found in 
ref. [H]. 

Finally we would like to emphasize that the qq Fock component of the photon represents the 
highest twist shadowing correction |T7], and vanishes at large quark masses as l/rrij. This does 
not happen for higher Fock states containing gluons, which lead to GS. Therefore GS represents 
the leading twist shadowing correction pU |44] . Moreover, a steep energy dependence of the 
dipole cross section aqq{r,s) (see Eq. ((61)) especially at smaller dipole sizes r causes a steep 
energy rise of both corrections. 



3 Gluon shadowing 

In the LC Green function approach [191 [TTl [20l [2T|, [18] the physical photon I7*) is decomposed 
into different Fock states, namely, the bare photon |7*)o, plus \qq), \qqG), etc. As we mentioned 
above the higher Fock states containing gluons describe the energy dependence of the photoab- 
sorption cross section on a nucleon, and also lead to GS in the nuclear case. However, these 
fluctuations are heavier and have a shorter coherence time (lifetime) than the lowest \qq) state, 
and therefore at small and medium energies only the \qq) fluctuations of the photon matter. 
Consequently, GS, which is related to the higher Fock states, will dominate at higher energies, 
i.e. at small values of xbj < 0.01. Since we will study the shadowing region of xsj ^0.01 and 
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the available experimental data reach values of xbj down to ~ 10~^, we will include GS in our 
calculations and show that it is not a negligible effect. Besides, no data for gluon shadowing 
are available and one has to rely on calculations. 

In the previous Section [2] we discussed the nuclear shadowing for the \qq) Fock component of 
the photon. It is dominated by the transverse photon polarizations, because the corresponding 
photoabsorption cross section is scanned at larger dipole sizes than for the longitudinal photon 
polarization. The transverse qq separation is controlled by the distribution amplitude Eq. ([8]), 
with the mean value given by Eq. (fT2!) . Contributions of large size dipoles come from the 
asymmetric qq fluctuations of the virtual photon, when the quark and antiquark in the photon 
carry a very large (a ^ 1) and a very small fraction [a 0) of the photon momentum, and 
vice versa. The LC wave function for longitudinal photons (1281) contains a term (1 — a)^, 
which makes considerably smaller the contribution from asymmetric qq configurations than 
for transversal photons (see Eq. (I27j) ). Consequently, in contrast to transverse photons, all qq 
dipoles from longitudinal photons have a size squared (x and the double-scattering term 
vanishes as oc l/Q"^. The leading-twist contribution for the shadowing of longitudinal photons 
arises from the \qqG) Fock component of the photon because the gluon can propagate relatively 
far from the qq pair, although the q-q separation is of the order After radiation of the 

gluon the pair is in an octet state, and consequently the \qqG) state represents a GG dipole. 
Then the corresponding correction to the longitudinal cross section is just gluon shadowing. 

The phenomenon of GS, just as for the case of nuclear shadowing discussed in the Introduc- 
tion, can be treated differently depending on the reference frame. In the infinite momentum 
frame this phenomenon looks similar to gluon-gluon fusion, corresponding to a nonlinear term 
in the evolution equation [15]. This effect should lead to a suppression of the small-xsj gluons 
also in a nucleon, and lead to a precocious onset of the saturation effects for heavy nuclei. 
Within a parton model interpretation, in the infinite momentum frame of the nucleus the gluon 
clouds of nucleons which have the same impact parameter overlap at small xbj in the longitu- 
dinal direction. This allows gluons originated from different nucleons to fuse, leading to a gluon 
density which is not proportional to the density of nucleons any more. This is gluon shadowing. 

The same phenomenon looks quite different in the rest frame of the nucleus. It corresponds 
to the process of gluon radiation and shadowing corrections, related to multiple interactions 
of the radiated gluons in the nuclear medium [28]. This is a coherence phenomenon known 
as the Landau-Pomeranchuk effect, namely the suppression of bremsstrahlung by interference 
of radiation from different scattering centers, demanding a sufficiently long coherence time of 
radiation, a condition equivalent to a small Bjorken xbj in the parton model. 

Although these two different interpretations are not Lorentz invariant, they represent the 
same phenomenon, related to the Lorentz invariant Reggeon graphs. It was already discussed 
in detail in refs. [201 SS] that the double-scattering correction to the cross section of gluon 
radiation can be expressed in Regge theory via the triple-Pomeron diagram. It is interpreted 
as a fusion of two Pomerons originated from different nucleons, 2IP — > IP, which leads to a 
reduction of the nuclear gluon density Ga- 

Notice that in the hadronic representation such a suppression of the parton density cor- 
responds to Gribov's inelastic shadowing [25], which is related to the single diffraction cross 
section. In particular, GS corresponds to the triple-Pomeron term in the diffractive dissociation 
cross section, which enters the calculations of inelastic corrections. 

There are still very few numerical evaluations of gluon shadowing in the literature, all of 
them done in the rest frame of the nucleus, using the idea from ref. [28]. As was discussed 
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above gluon shadowing can be identified as the shadowing correction to the longitudinal cross 
section coming from the GG dipole representing the \qqG) Fock component of the photon. 
An important point for the evaluation of GS is knowing about the transverse size of this GG 
dipole. This size has been extracted in ref. [24] from data for diffractive excitation of the incident 
hadrons to the states of large mass, the so called triple-Pomeron region. The corresponding 
diffraction cross section (cx r^) is a more sensitive probe of the mean transverse separation than 
the total cross section (oc r^). Consequently, it was found in ref. [23] that the mean dipole 
size of the GG system (radius of propagation of the LC gluons) is rather small , tq ~ 0.3 fm 
[47] . Such a small quark-gluon fluctuation represents the only known way how to resolve the 
long-standing problem of the small size of the triple-Pomeron coupling. 

To incorporate the smallness of the size of quark-gluon fluctuations into the LC dipole ap- 
proach, a nonperturbative LC potential describing the quark-gluon interaction was introduced 
into the Schrodinger equation for the LC Green function describing the propagation of a quark- 
gluon system. The strength of the potential was fixed by data on high mass {M\) diffraction 
pp — > pX [21]. This approach allows to extend the methods of pQCD to the region of small 
Q^. Since a new semihard scale 1/ro ~ 0.65 GeV is introduced, one should not expect a sub- 
stantial variation of gluon shadowing at < 4/rQ. Indeed, the calculations performed in [23] 
for = and 4 GeV^, using different techniques, led to about the same gluon shadowing. At 
higher shadowing slowly (logarithmically) decreases, in accordance with the expectations 
based on the evolution equation [3], which clearly demonstrates that GS is a leading-twist effect. 

In this paper we repeated the calculations [2^ of the ratio of the gluon densities in nuclei 
and nucleon, 

, 2._ Ga{xbj,Q^) _ AatotiqqG) 

^ ) - AG,ixs„Q^) - ' - ' ^''^ 

where AatotiqqG) is the inelastic correction to the total cross section atot"^, related to the 
creation of a \qqG) intermediate Fock state. 



AatotiqqG) = Re J dz2 J dzi pAib, zi) pAib, Z2) J d'^X2d'^y2d'^xid'^yi J da. 



d ac 

q 



X F^,_,ggQix2, m, oiq, oig) GggGix2, m, Z2] fi, ^i, Zi) F^^^ggdxi, yi, ag, ac) ■ (56) 

Here x and y are the transverse distances from the gluon to the quark and antiquark, respec- 
tively, ag is the fraction of the LC momentum of the qq carried by the quark, and ac is the 
fraction of the photon momentum carried by the gluon. F^*^gqG is the amplitude of diffractive 
qqG production in a 7*A^ interaction [23], and it is given by 
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F^,^ggGix,y,aq,aG) = -'^ ggiag,x - y) 



X 



(yqqix) + (Tqqiy) - ^^^(f - y) 



(57) 



where '^qq and '^gc are the LC distribution functions of the qq fluctuations of a photon and qG 
fluctuations of a quark, respectively. 

In the above equation GqqGix2,y2, Z2; Xi,yi, Zi) is the LC Green function which describes 
the propagation of the qqG system from the initial state with longitudinal and transverse 
coordinates zi and xi,yi, respectively, to the final coordinates (-22, a^2, ^2)- For the calculation 
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of gluon shadowing one should suppress the intrinsic qq separation, i.e. assume x = y. In this 
case the Green function simphfies, and effectively describes the propagation of a gluon-gluon 
dipole through a medium. 

An important finding of ref. [24j is the presence of a strong nonperturbative interaction 
which squeezes the gluon-gluon wave packet and substantially diminishes gluon shadowing. The 
smallness of the gluon-gluon transverse separation is not a model assumption, but is dictated 
by data for hadronic diffraction into large masses (triple- Pomeron regime), which is controlled 
by diffractive gluon radiation. 
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Figure 2: The ratio of the nucleus-to-nucleon gluon densities as function of 
the thickness of the nucleus, L = T{b)/po, at = 4 GeV^ and different fixed 
values of xbj- Figure is taken from ref. [20] . 



Further calculational details can be found in [21]. In our case we calculated the gluon 
shadowing only for the lowest Fock component containing just one LC gluon. In terms of 
the parton model it reproduces the effects of fusion of many gluons to one gluon (in terms of 
Regge approach it corresponds to the nTP IP vertex). Inclusion of higher multigluon Fock 
components is still a challenge. However, their effect can be essentially taken into account 
by the eikonalization of the calculated RdxBjyQ'^), as argued in [IS]. In other words, the 
dipole cross section, which is proportional to the gluon density at small separations, should be 
renormalized everywhere, in the form 

(^qq =^ Rg (^qq ■ (58) 

Such a procedure makes the nuclear medium more transparent. This could be expected since 
Gribov's inelastic shadowing is known to suppress the total hadron-nucleus cross sections, i.e. 
to make nuclei more transparent [29| HQ] . 
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As an illustration of not very strong onset of GS, here we present Rc{xBj,Q'^), Eq- dSS]), 
for different nuclear thicknesses TA{b). Using an approximation of constant nuclear density 
(see Eq. ( I36i) ). Ta(6) = po L, where L = 2 \J R\ — 6^, the ratio Rcixsj, Q"^) is also imphcitly a 
function of L. An example for the calculated L-dependence of Rcixsj, Q"^) at = 4 GeV^ is 
depicted in Fig. [2] for different values of Xsj- 

One can expect intuitively from Eq. (1581) that GS should always diminish the nuclear cross 
sections of various processes in nuclear targets, and that the onset of GS is stronger for heavier 
nuclei. However, this is not so for incoherent electroproduction of vector mesons, analyzed in 
ref. [20] . The specific structure of the expression for the nuclear production cross section causes 
that the cross section of incoherent electroproduction of vector mesons is rather insensitive 
to GS. Furthermore, the effect of GS is stronger for light than for heavy nuclear targets, in 
contradiction with the standard intuition. Moreover, for heavy nuclei the effect GS can lead 
even to a counterintuitive enhancement (antishadowing), as was analyzed in ref. [20]. For the 
case of coherent vector meson production 7*^4 — >■ V A [20], GS was shown to be a much stronger 
effect in comparison with incoherent production, which confirms the expected reduction of the 
nuclear production cross section. 

Similarly, it was analyzed in ref. [16] that multiple scattering of higher Fock states containing 
gluons leads to an additional suppression of the Drell-Yan cross section. In the present paper we 
will demonstrate that gluon shadowing also suppresses the total photoabsorption cross section 
on a nucleus a'tot^{xBj, Q^)- Here we expect quite a strong effect of GS in the shadowing region 
of small XBj ^ (0.01 -v- 0.001), in the kinematic range of available data corresponding to small 
and medium values of ~ a few GeV^. 

4 Numerical results 

As we mentioned above the main goal of this paper is to compare for the first time available 
experimental data with realistic predictions for nuclear shadowing in DIS, based on exact 
numerical solutions of the evolution equation for the Green function. Such a comparison is 
performed for the shadowing region of small xbj ^ 0.01. As was discussed in the previous 
section one should take into account also a contribution of gluon shadowing, which increases the 
overall nuclear suppression. The effect of GS was already calculated in ref. [50j, but only for the 
F2 / F2 ratio of structure functions. Although the quark shadowing was compute approximately 
via longitudinal form factor of the nucleus and assuming only the leading shadowing term, it 
was shown that GS is a rather large effect at xsj ~ 10~^. In the present paper we will also 
show that GS is not a negligible effect, and can in principle be detected by the data on the 
total photoabsorption nuclear cross section in future experiments. 

The predictions for nuclear shadowing based on an exact numerical solution of the evolution 
equation was compared in ref. [18] with approximate results obtained using the harmonic os- 
cillatory form of the Green function (fT6l) . Quite a large discrepancy was found, in the range of 
XBj ^ 0.001, where the variation of the transverse size of the qq pair during propagation through 
the nucleus becomes important. Such a variation is naturally included in the Green function 
formalism and consequently the exact shape of the Green function is extremely important. 

We use an algorithm for the numerical solution of the Schrodinger equation for the Green 
function, as developed and described in ref. [18] . This gives the possibility of calculating nuclear 
shadowing for arbitrary LC potentials V^g(z,r, a) and nuclear density functions. Because the 



17 



available data from the E665 [51], [521 and NMC [531 ISl] collaborations cover the region of 
small and medium values of < 4 GeV^, we prefer the KST parametrization of the dipole 
cross section ([6]), which is valid down to the limit of real photoproduction. On the contrary, 
the second GBW parametrization [32] of the dipole cross section cannot be applied in the 
nonperturbative region and therefore we do not use it in our calculations. 

In the process of exact numerical solutions of the evolution equation for the Green function, 
the imaginary part of the LC potential contains the corresponding KST dipole cross section 
as well. The nuclear density function pA{b, z) was taken in the realistic Wood-Saxon form, with 
parameters taken from ref. [43]. The nonperturbative interaction effects between q and q are 
included explicitly via the real part of the LC potential of the form ( JTSll . which is supported 
also by the fact that the data from E665 and NMC collaborations correspond to very small 
values of < 1 GeV^ in the region of small xsj < 0.004. 

We included also the effects of gluon shadowing for the lowest Fock component containing 
just one LC gluon. Although the inclusion of higher Fock components with more gluons is 
complicated, their effect was essentially taken into account by eikonalization of the calculated 
RcixBjyQ"^) [1H|, i-e using the renormalization ( l58l) . 

Nuclear shadowing effects were studied via the xbj- behavior of the ratio of proton structure 
functions (H3l) divided by the mass number A. First we present nuclear shadowing for a lead 
target in Fig. [3] at different fixed values of Q^. The thick and thin solid curves represent the 
predictions obtained with and without the contribution of gluon shadowing, respectively. 

One can see that the onset of GS happens at smaller xbj than the quark shadowing, which 
is supported by the fact that higher Fock fluctuations containing gluons are in general heavier 
than qq and have a shorter coherence length. Fig. [3] demonstrates quite a strong effect of GS 
in the range of xbj G (0.01,0.0001) where the most of available data exist. This is a result of 
the suppression of the dipole cross section by the renormalization fl58l) . which can result only 
in a reduction of the total photoabsorption cross section on a nuclear target. Besides the effect 
of GS is stronger at smaller because corresponding Fock fluctuations of the photon have a 
larger transverse size. 

In Fig.[3]we also present, for comparison and by the dotted lines, the approximate predictions 
for nuclear shadowing in DIS using constant nuclear density ( !36l) and the quadratic form of 
the dipole cross section, o"(r, s) = C{s) r^. The energy dependent factor C(s) is determined by 
Eq. (jllj), and the uniform nuclear density is fixed by the condition ( H5ll . One can see that these 
approximate predictions overestimate the values of nuclear shadowing obtained by means of 
an exact numerical solution of the evolution equation for the Green function. The difference 
from the exact calculation (thin solid lines) is not large and rises towards small values of Q^. 
The reason is that the quadratic approximation of the dipole cross section cannot be applied 
exactly at large dipole sizes. Since the available data from the E665 [511 112] and NMC [521 IM] 
collaborations at smallest values of xbj correspond also to small -C 1 GeV^, one can expect 
a larger difference between the exact and approximate results in comparison with what is shown 
in Fig. [3] at = 2 GeV^. Keeping the quadratic form of the dipole cross section, but using 
the realistic nuclear density [13] , one can obtain the results depicted in Fig. [3] by the dashed 
lines. It brings a better agreement with the exact calculations. 

At low XBj ^ 10""^ one should expect a saturation of nuclear shadowing at the level given 
by Eq. (HTj) . This is realized only for the dipole cross section, without energy dependence, i.e. 
for example for parametrization fl35|) of the dipole cross section with constant factor C{s) ^3 
[To] . However, this is not so for the realistic KST parametrization Eq. ([6]), where the saturation 
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level is not fixed exactly due to energy (Bjorken xbj-) dependence of the dipole cross section 




-1 



Figure 3: Nuclear shadowing for lead. Calculations correspond to exact 
numerical solution of the evolution equation for the Green function using 
the KST [21] parametrization of the dipole cross section and a realistic 
nuclear density function of the Woods-Saxon form ^Sj. The thick and 
thin solid curves represent the predictions calculated with and without 
contribution of gluon shadowing, respectively. The dotted lines are cal- 
culated using a constant nuclear density function fl5Bl) and the quadratic 
form of the dipole cross section, cr{r,s) = C(s)r^, where the energy de- 
pendent factor C(s) is determined by Eq. fHi|) . The dashed curves are 
calculated for the same quadratic form of the dipole cross section, but 
for the realistic nuclear density function ([43j) 
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Figure 4: Comparison of the model with experimental data from the 
E665 [HH [52] and NMC [53], [5^ collaborations. Calculations correspond 
to the exact numerical solution of the evolution equation for the Green 
function using KST [24J parametrization of the dipole cross section and 
a realistic nuclear density function of the Woods-Saxon form [43]. The 
solid and dashed curves are calculated with and without the contribution 
of gluon shadowing, respectively. 



In Fig. m we present a comparison of the model predictions with experimental data at 
small Xfij, from the E665 (SH [S3 and NMC [531 [M] collaborations. One can see a quite 
reasonable agreement with experimental data, in spite of the absence of any free parameters in 
the model. Several comments are in order: first, if OS is not taken into account, for the C/D 
and Ca/D ratios the nuclear shadowing looks overestimated in comparison with the E665 data 
for XBj ~ 0.01, while it looks in a good agreement for C/D, and a little bit underestimated 
for Ca/D in comparison with the NMC data. This is affected by the known incompatibility 
of the results from both experiments for the ratios over D. Second, as was discussed in the 
previous section, the effect of GS produces an additional nuclear shadowing which rises with 
mass number A. Consequently, it leads to a small overestimation of the nuclear shadowing for 
the C/D and Ca/D ratios in comparison with the E665 data, but it seems to be a in good 
agreement with the NMC data. 
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For heavy nuclear targets there is only E665 data for the ratios Xe/D and Pb/D. Fig. S] 
shows a reasonable good description of these data, even if the effect of GS is taken into account. 
The difference between the solid and dashed lines in Fig. H] represents quite a large effect of 
GS, which was neglected up to the present time in calculations of nuclear shadowing in DIS 
[T9t \T7\ [23] assuming that it would be a very small effect in the kinematic range covered by the 
available experimental data. On the contrary, looking at Fig. Hlone can see that the effect of GS 
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Figure 5: Model predictions for nuclear shadowing for a broad x^j-region down to 
10~^ corresponding to LHC kinematical range at three different values of = 2, Q 
and 18 GeV^. Calculations of the nuclear shadowing for the qq Fock component of the 
photon correspond to the exact numerical solution of the evolution equation for the Green 
function using KST [21] parametrization of the dipole cross section and a realistic nuclear 
density function of the Woods-Saxon form j33]. The solid and dashed curves are calculated 
with and without the contribution of gluon shadowing, respectively. 
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Figure 6: The same as Figl5] but with GBW [32] parametrization of the dipole cross 
section. 



as an additional nuclear shadowing cannot be neglected and should be included in calculations 
already in the region of xbj < 0.01 0.001. Very large error bars especially at small xbj ~ 10~^ 
do not allow to investigate separately the effect of GS, and therefore more exact new data on 
nuclear shadowing in DIS at small xbj are very important for further exploratory studies of the 
nuclear modification of structure functions and also for gluon shadowing. 

For completeness we present also in Figs. E] and [H] predictions for nuclear shadowing down to 
very small xbj = 10"'' accesible by experiments at LHC using two different realistic parametriza- 
tions of the dipole cross section, KST [2^ and GBW [32] . Again, one can see quite large effect 
of GS as a difference between the solid and dashed lines. 

Here, we would like to emphasize that at xbj ^ 10~^ transverse size separations of the 
photon fluctuations are "frozen" during propagation through nuclear medium and one can use 
the simplified expressions, Eq. fl40l) and fl41l) for calculation of nuclear shadowing. 
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Figure 7: Comparison of the model results for the ratio Pb/nucleon 
obtained without (upper thick solid line) and with (lower thick solid line) 
gluon shadowing, with other models, versus XBj, at fixed = 3 GeV^. 
Bartels are the results from [58J, Frankfurt from [SS] {Q^ = 4 GeV^), and 
ACKLS from pj. 



Finally, we present in Fig. [7] a comparison of the nuclear shadowing calculated using our 
model with the results of other models, for = 3 GeV^ (except the results of ref. which 
are at = 4 GeV^). Notice that the difference between models rises towards small values of 
xqj^ as a result of the different treatment of various nuclear effects, and absence of relevant 
experimental information at such small xsj- At xsj = 10^^ we predict quite a large effect of 
GS (compare upper and lower thick solid lines). 

In ref. [56] nuclear structure functions were studied using relation with diffraction on nu- 
cleons known as Gribov inelastic corrections. The results of these calculations are depicted in 
Fig. El by dotted curves. 

The model presented in |55] employs again a parametrization of hard diffraction at the scale 
Qq, which gives nuclear shadowing in terms of Gribov's corrections similar to ref. |56j- Then 
the nuclear suppression calculated at Qq is used as initial condition for Dokshitzer-Gribov- 
Lipatov-Altareli-Parisi (DGLAP) [57J evolution. This results are presented in Fig. [7] by dashed 
curves. 
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The model based on a numerical solution of the non-linear equation for small-x^j evolution 
in nuclei was employed in [5^. The result is shown in Fig. [7] by thin solid curve. 

At very small xbj our model predictions, including the effect of GS, roughly agrees with 
those of [56], but lie below the results of other models. If the effect of GS is not taken into 
account the situation is substantially different and the corresponding curve (see upper thick 
solid line) lies in between the results from other models. 

Most models presented above are based on eikonal formulas, which should be used only in the 
high energy limit, when the coherence length ^ Ra, at xbj ^ 10~^. However, they were 
applied also in the region when Ic < Ra- In this transition shadowing region xbj G (0.0001, 0.01) 
such approximations lead in general to a larger nuclear shadowing than a realistic situation, 
when more exact expressions should be more appropriate (compare Eqs. (l39l) and ( HTl) ). For 
this reason, theoretical predictions of most models overestimate nuclear shadowing in the range 
of XBj where available experimental data exist. 

So far the main source of experimental information on gluon shadowing was DIS on nuclei. 
Although it probes only quark distributions, the dependence of nuclear effects is related via 
the evolution equations to the gluon distribution. At the smallest value of xbj = 0.01 reached in 
the NMC experiment the gluon suppression factor Rc^Sn) / Rq{C) = 0.87 ± 0.05 was obtained 
in ref. [59] within the Leading-Log (LL) approximation. This result is somewhat lower that 
our expectation Rq ~ 0.98 which can be read out from Fig. [2l However, according to [59] 
the next-to-LL corrections at xbj = 0.01 are about 10 — 20%, which apparently eliminates the 
disagreement with our calculations. Furthermore, the full Leading-Order (LO) DGLAP analysis 
of the NMC data [53', 'M] in ref. [60], which should not be less accurate than LL calculations, led 
to a conclusion that the NMC data are not sensitive to gluon shadowing. Moreover, the recent 
Next-to-LO (NLO) analysis by de Florian and Sassot [61| was claimed to be sensitive to gluons. 
This analysis found almost no GS at xbj = 0.01 in good agreement with our calculations. 

Other possible sources of information about gluon shadowing were considered in refs. [S21 
ESI El] • It was proposed in ^2] to probe gluons in nuclei by direct photons produced in p — A 
collisions in the proton fragmentation region where one can access smallest values of the light- 
front momentum fraction variable X2 in nuclei. This, however, should not work, since at large 
value of the light-front momentum fraction variable xi in the proton (i.e. at large Feynman 
xf) one faces the energy sharing problem [65]: it is more difficult to give the whole energy 
to one particle in p — A, than in p — p collision. This effect leads to a breakdown of QCD 
factorization and to nuclear suppression observed at forward rapidities [651 EH [671 (Mj in any 
reaction measured so far, even at low energies, where no shadowing is possible. 

An attempt to impose a restriction on GS analyzing the nuclear effects in J/\I' production 
observed in p — A collisions by the E866 experiment [HH] , was made in [03] . The mechanisms 
of J/"^ production and nuclear effects are so complicated, that it would be risky to rely on 
oversimplified models. Indeed, the analysis performed in this paper completely misses the color 
transparency effects, which are rather strong [70] and vary throughout the interval ofxp studied 
in this paper. For this reason the results of the analysis are not trustable. 

The new analysis of nuclear parton distribution functions performed in [M] included the 
BRAHMS data for high-p^- pion production at forward rapidities |7Ij. As we mentioned above, 
hadron production in this kinematic region of large xi (xp) is suppressed by multiple parton 
interactions |65j, rather than by shadowing. Consequently, the results of this analysis are 
not trustable either. Moreover, it seems to provide another confirmation for an alternative 
dynamics for the suppression observed in the data [71]. Indeed, it was concluded in [61] that 
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gluons in lead target are completely terminated at XBj = 10~^ where Rq < 0.05 is predicted. 
This is cannot be true because in the limit of strong shadowing the gluon ratio has a simple 
form Rg = TrR\/{Aaeff), where a^ff is the effective cross section responsible for shadowing. 
The strong effect predicted in [64] needs cTe// > 150 mb. 

5 Summary and conclusions 

We presented a rigorous quantum-mechanical approach based on the light-cone QCD Green 
function formalism which naturally incorporates the interference effects of CT and CL. Within 
this approach HH [TTl [IB] we studied nuclear shadowing in deep-inelastic scattering at small 
Bjorken xbj- 

Calculations of nuclear shadowing corresponding to the qq component of the virtual photon 
performed so far were based only on efforts to solve the evolution equation for the Green function 
analytically, and unfortunately an analytical harmonic oscillatory form of the Green function 
( fT6l) could be obtained only by using additional approximations, like a constant nuclear density 
function ( |36ll and the dipole cross section of the quadratic form (1351) . This brings additional 
theoretical uncertainties in the predictions for nuclear shadowing. In order to remove these 
uncertainties we solve the evolution equation for the Green function numerically, which does 
not require additional approximations. 

In ref. [H] it was found for the first time the exact numerical solution of the evolution 
equation for the Green function, using two realistic parametrizations of the dipole cross section 
(GBW [32] and KST [S^), and a realistic nuclear density function of the Woods-Saxon form [13]. 
It was demonstrated that the corresponding nuclear shadowing shows quite large differences 
from approximate results [T9l [T7] . On the other hand, we showed that approximate calculations 
corresponding to uniform nuclear density (HSi) and quadratic dipole cross section (H7|) , but with 
energy dependent factor C{s) determined by Eq. fl44l) . bring a better agreement with exact 
realistic calculations (see Fig. [3]). However, the difference from the exact calculations rises 
towards small values of Q^, where available data exist at smallest values of xbj ~ 10~^. This 
confirms the claim that the quadratic approximation of the dipole cross section cannot be 
applied at large dipole sizes. 

Since the available data from the shadowing region of xbj ^ 0.01 comes mostly from the 
E665 and NMC collaborations, and cover only small and medium values oi <A GeV^, we 
used only the KST realistic parametrization [23] of the dipole cross section, which is more 
suitable for this kinematic region and the corresponding expressions can be applied down to 
the limit of real photoproduction. On the other hand, the data obtained at the lower part 
of the x^j-kinematic interval correspond to very low values of < 1 GeV^ (nonperturbative 
region). For this reason we include explicitly the nonperturbative interaction effects between 
q and q, taking into account the real part of the LC potential Vqq f|T5|) in the time-dependent 
two-dimensional Schrodinger equation f|T3l) . 

In order to compare the realistic calculations with data on nuclear shadowing, the effects of 
GS are taken into account. The same path integral technique [23] can be applied in this case, and 
GS was calculated only for the lowest Fock component containing just one LC gluon. Although 
the inclusion of higher Fock components containing more gluons is still a challenge, their effect 
was essentially taken into account by eikonalization of the calculated RcixBjjQ"^) ,using the 
renormalization (l58l) . We found quite a large effect of GS, which starts to be important already 
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at XBj ~ 0.01. The effect of GS rises towards small xbj because higher Fock components with 
more gluons having shorter coherence time will contribute to overall nuclear shadowing. Such 
a situation is illustrated in Fig. [31 

Performing numerical calculations, we find that our model is in reasonable agreement with 
existing experimental data (see Fig.H]). Large error bars and incompatibility of the experimental 
results from the E665 and NMC collaborations do not allow to study separately the effect of 
GS, and therefore more accurate new data on nuclear shadowing in DIS off nuclei at still smaller 
XBj ^ 10~^ are very important for further exploratory studies of GS effects. 

Comparison among various models shows large differences for the Pb/nucleon ratio of struc- 
ture functions at xbj = 10~^ and = 3 GeV^ (see Fig. [Tj), which has a large impact on the 
calculation of high-p^ particles in nuclear collisions at RHIC and LHC. Such large differences at 
small XBj among different models should be testable by the new more precise data on nuclear 
structure functions, which can be obtained in lepton-ion collider planned at BNL |72j . 

In most models presented above the final formulae for nuclear shadowing are based on the 
eikonal approximation, which can be used exactly only in the high energy limit, Ic ^ Ra- 
Consequently, such an approach cannot be really applied in the transition shadowing region, 
Ic ~ Ra, where XBj G (0.0001,0.01), because it produces a larger nuclear shadowing than 
in a realistic case when more appropriate expressions should be taken into account (compare 
Eqs. (139]) and dUD). 

Concluding, a combination of the exact numerical solution of the evolution equation for the 
Green function with the universality of the LC dipole approach based on the Green function 
formalism provides us with a very powerful tool for realistic calculations of many processes. 
It allows to minimize theoretical uncertainties in the predictions of nuclear shadowing in DIS 
off nuclei, which gives the possibility to obtain reliable information about nuclear modification 
of the structure functions at low Xsj, with an important impact on the physics performed in 
heavy ion collisions at RHIC and in lepton-ion interactions planned at BNL. 
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